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A THEORETICAL ANALYSIS OF THE 
MSFC HOLOGRAPHIC CORRELATION SYSTEM 

Abstract 

A theoretical analysis has been completed for the correlation 
output signal for the MSFC holographic correlation filtering system. 
Under appropriate assumptions, the correlation output was derived as a 
function of the roughness of the tested surface, the displacement of 
the illuminated area of the test surface, the characteristics of the 
optical components used in the system, and the system configuration. 

In addition, an approximate relationship between the displacement of the 
detected signal (which is focused on the photomultiplier tube) 

and the displacement of the illuminated area on the test surface was 
also derived. 


A THEORETICAL ANALYSIS OF THE 
MSFC HOLOGRAPHIC CORRELATION SYSTEM 

I. Introduction 

Holographic non-destructive testing (HNDT) has a broad range of 
applications in the military and industrial testing of ultra-expensive 
and zero-defect- demanded instruments and equipment such as the main 
shuttle engine of NASA's space shuttles. 

Most of the work of HNDT has been done in holographic inter- 
ferometry, which basically involves recording and interpreting fringes 
resulting from the interaction between two mutually coherent wavefronts. 
Recently, a practical method of interpreting the interference fringes 
has been devised end experimentally verified . 1 * 2 The method enables one 
to determine anomalous displacement on the surface of an object semi- 
quantitatively. On the other hand, the technique of holographic corre- 
lation can be used to zero-in on a specific abnormal region, which can 
be determined by the previous holographic interferometry method. The 
correlation technique will give an overall number which signifies the 
similarity of the two wavefronts being processed. From this number, the 
knowledge of the stress and strain on that specific region can be found. 

The correlation technique utilized a Vander Lugt filter as a 
matched filter in the study of the autocorrelation properties of the 
stressed and unstressed object . 3 * 8 The main reason which motivated the 
contracted study described in this report is due to the experimental 
discovery a secondary peak in the putput of the correlation intensity 
in the MSFC corelation system. In addition, it was found that there 
was no variations in the correlation signal strength when an aperture 
was inserted in the object beam. The anomalism of the secondary peak 
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was removed in -the experiment after the contractor suggested that the 
system configuration be changed to the present one with an oblique 
incident object beam from its previous normal incident setup. Moreover, 
an object beam focusing lens was added to focus the beam on the rough 
test surface. The focused beam causes the scattered light from the sur- 
face to have a broader speckle pattern. A hologram made in the present 
system was able to produce a correlation function with a much stronger 
signal-to-noise ratio (about two orders of magnitude stronger than the 
signal-to-noise ratio obtained in the previous system). Results showed 
that no more secondary peak was observed. Hence the secondary peak 
observed before probably was due to the interference of noise to an 
extremely weak correlation signal. The effect of aperture, however, 
should be easily visualized after the theoretical analysis is presented 
in the following. 

Section II will be the statement of the problem and the basic 
assumptions being made. The theoretical analysis will be provided in 
Section III. Section IV are the conclusions. 

II. Statement of the problem and basic assumptions 

A. State of the problem 

A theoretical analysis is required to describe the correlation 
signal, detected in the holographic correlation filtering system as 
shown in Figure 1 , as a function of the roughness of the test surface, 
the displacement of the illuminated area, the characteristics of the 
optical components, and the configuration of the system. 

B. Basic assumptions 

It. is assumed that (1) The test surface is uniformly and iso- 
tropically rough in two dimensions with normally distributed roughness 
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D iffusive surface 


Point on 
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Figure 1. A diagram showing the wave vectors of the 

incident and reflected beams and the diffusive 
surface. 
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C (x*y) of zero mean ("flat")* variance o 2 , and correlation distance X. 
(2) Both a and X are of similar magnitude as the optical wavelength X 
or larger. (3) The radius of the laser beam* r 0 , focused on the object 
surface is large compared with the correlation distance X and the wave- 
length X. 


III. Theoretical Analysis 
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When the laser beam strikes a point P at the surface of the object, 
light will be scattered. As shown in Fig. 1, the scattering angles ©j , 

02 and 03 at P may be defined according to the Bechmann's convention. 9 

*•> 

Is the angle of incidence of the beam with wave vector k L on the 
mean surface, 

0i - cos" 1 [ - (n • ki) / |k! ( ], (1) 

^ -> 

where [ki| = 2 n/A which is the absolute value of ki and n is the 
unit vector of the outer normal through the point P of the mean surface. 
Similarly, 

e 2 * cos" 1 [ (n • k 2 ) / |k 2 | ], (2) 

•> 

where k 2 is the wave vector, of the scattered field. 

e 3 is the angle between the planes formed by n, k 2 and n, k 2 , hence 

e 3 = cos" 1 [ (n x ki) • (r. x k 2 ) / (s i ne j s1ne 2 ) ], (3) 

where 6i 1 0 and e 2 i 0. 

If P is at the origin, these angles are defined by 010 , 8 2 q, and 
030 respectively. 

In addition, a scattering vector can be defined by 

♦ * 2TT r / . * 

V ■ ki - k 2 = ~ [ ( s in 0i 0 - sine 20 cose 30 ) x 
- ( sine 20 sine 30 ) y - (cose 10 + cose 20 ) * ], (4) 

where x, y, z are unit vectors in the x, y, z directions, respectively. 

On the rough surface, the true local scattering angles are in 
general different from ©i, e 2 , e 3 . Let s be the true local normal unit 
vector which can be expressed by the slopes ; x * at/3x and cy * 8c/3y 
of the roughness c (x»y): 

*s = ( -tx*x “ Cy y + z ) (1 + cj + c 2 


(5) 
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The local scattering angles can be defined by replacing n by s in 
Equations (1) through (3). 

Following Bechman 1 , the scattered field at the far distance R 0 from 
the surface after adaption to a Gaussian beam can be written as 


Ao — 


exp [ 2l ^ 01 ] cos 9io 


2n 2 r 0 2 R 0 


/ / P (x,y) F (c x , c y ) • 


e^ v ' r d x d y, (6) 

where P (x,y) = exp [ - (x 2 cosQ 10 + y 2 ) / r 0 2 ], (7) 

and F (c x , Sy) is defined in terms of the reflection coefficient R and 
the scattering angles 0 l5 0 2 , and 0 3 : 

F U x , Cy) = [ (1 - R) s i n 0 ! + (1 + R) sin0 2 cosa 3 ] r, x 

+ [ (1 + R) sin0 2 sin0 3 ] t,y 

- (1 + R) cos 0 2 + (1 - R) coso ! . (8) 

Since on the rough surface, the true local polarization state of 
the incident field is in general unknown, F (c x , cy) can only be treated 
as a random variable and statistical analysis will be required. 


The optical phase v • r in Eq. (6) can be written as 
v • r = v x x + Vy y + v z (z + e) 

= <t> (x,y) + v z (z + c ), (9) 

where $ (x,y) = v x x + Vy y . (10) 

If a hologram is made in the configuration as shown in Fig. ? and 
replaced to its original position, the hologram serves as a matched filter 
to the object beam when the reference beam is turned off. If the correla- 
tion ] en c; i s not in the system, the effect of the amplitude transmitted 
by the hologram at a distance R 0 from P 0 is proportional to the follow- 
ing correlation function with the symbol < > denoting a statistical 
ensemble average. 
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From Figure 2. A diagram of the laser holographic 

Spectra Physics correlation filtering system. 

Model 125A 
HeNe laser 
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6l (av) »‘<A 2 (v) A 2 * (v + av) > , (It) 

Substituting Eqs. (6) and (9) into Eq. (11)* one obtains 

k 2 cos 2 ajo * 

— exp [ i (y zl Zi - v z2 z 2 ) ] ffff P ( x j * yj 

4n 4 r 0 - Ro 2 - 

P (x 2 , y 2 ) exp { 1 [ 0 (xi # yi) - ♦ (x 2 , y 2 ) ] ) 

X < F (c x i » ; y i) F (c x2 , C y2 ) exp [ i (v 21 Ci - v z2 t 2 ) ] > 

dx t dy j dx 2 <ty 2 , (12) 

In Eq. (12), c 2 » C x2 ® nd Cy 2 have to be evaluated at x 2 - D x , 
y 2 - fy, if 5i» C x i» and Cyi are taken at Xj, yj. To calculate the 

ensemble average as an expectation value, the joint probability density 

function of the six normal random variables c 2 , c x i * c x2 » Cyi ♦ c y2 

is required. Because ;, ; x , and Cy are statistically independent, the 

density factorizes into three joint densities for the variable pairs 

(ci. C 2 )» (cxi. C x2 ) > Uyi* Cy 2 ). Each density is determined by the 

second moments of the variables. 10 

If the second moment of the roughness is given by 

< Ci C 2 > ■ O 2 C (x lt y lt x 2 - D x , y 2 - Dy) , (13) 

with C the autocorrelation coefficient of the roughness. 

Through some mathematical maniputation^one can write (See Appendix A), 
*> ■* 

Gj (A v, D) * n 2 exp [ - (o A v z ) 2 / 2 ] 
x /// P (xj, yi) P (x 2 , y 2 ) 

x exp { - q 2 [ ( x j - x 2 + D x ) 2 + (yi - y 2 + D y ) 2 ] 

+ 1 C ♦ (xj, yi) - $ (x 2 , y 2 )] } d X j d yl d x2 dy 2 . (H) 

Where in the above equation, 
o ■ k cos®io ( < F 2 > )** 


/ ( 2 n 2 r 0 2 Ro), 


(15) 


< F 2t F 2 > ■ < F (c x i. Cyi) F (c X 2» ?y2> > 
■ < F 2 > 

* 4 < R 2 > (I ♦ COS0 1O COS0 2 O 
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- s1ne 10 s1na 2 o coso 3 o) 2 / 

( COS 0 j o + coso 20 ) 2 , (16) 

and < R 2 > is the mean square of the reflection coefficient. 

After the integration in Eq. (14) has been carried out* the result 

is a product of many exponential factors. Among these factors, only 

those terms that can noticeably deviate from unity are considered as 

* 

significant and retained. Consequently, 

Gi (a v, D) * [ n 2 r 0 2 n 2 / (*T S hi h 2 cose 10 )3 


exp [ - i (v x D x + v y D y ) ] . (17) 

with 

S ■ exp { - [ vj / 2h j 2 ♦ v| / 2 h 2 2 

♦ 1/8 r 0 2 (a v x 2 / cos 2 e 10 + a v y 2 ) 

+ 1/2 (o a v 2 ) 2 ] } , (18) 

T ■ exp [ (- D x 2 cos 2 e 10 - D y 2 ) / (2 r 0 2 ) ] , (19) 

hi 2 - 2q 2 + cos 2 e 10 / r 0 2 , (20) 

h 2 2 - 2q 2 + 1 / r 0 2 , (21) 


and v x * H (v xl + v x2 ) . a v x « v x2 - v xl , etc. The term q 2 is described in 
Appendix A. 

Furthermore, the correlation function for the intensity may be 
written as 

G 2 (a v, D) • [n 2 r 0 2 n 2 / (hi h 2 cosei 0 ) 3 

[ $ 2 (0) + S 2 (a v) T 2 3 . (22) 


where 


S (0) * exp [ - V 2 / (2h j 2 ) - / (2h 2 2 ) 3 . 


(23) 



For a special case where & v * 0, Eq. (22) becomes 

G 2 (0, D) - [ n 2 r 0 2 n 2 / (h 2 h 2 coso 10 ) ] S 2 (0) [ 1 + V ] 

• K 2 { 1 + exp [ — (D x 2 cos 2 e 10 + D y z ) ] } (24) 

2r 0 2 

with 

K 2 *'[ n 2 r 0 2 fl 2 / (hj h 2 coso 10 ) ] S 2 (U), (25) 

From Eq. (24) » It can be seen that Gg^.D) Is a function of the angles of 

the Incident beam and the reflected beam, the radius of the beam on the 

rough surface, the roughness of the surface, and the displacements D x , 
and D . When D z varies, the radius of the beam rg will vary accordingly. 
Hence the variation of D z will affect the magnitude of the correlation 
function as well. , 

After the correlation lens Is inserted, the scattered light from the 
hologram is focused at the detector surface inside the photo-multiplier 
tube. The Fourier transformation performed by, the correlation lens 
undoes the Fourier transformation by the transform lens to the signal, 
therefore, Equations (17) and (24) can be used to describe the correla- 
tion intensity output. 

When the object Is translated from its original position as shown 
by the dotted line of Figure 3, another phenomenon will occur. The 
focal point at the photo-multiplier tube will also change its position 
due to the refraction of the laser light by the lenses the diffrac- 
tion of the beam by thr hologram and the lenses. The original path of 
the laser beam and the path of the light after the displacement is 
shown by the solid and dashed lines respectively. For example, if 
the object is translated along the Incident laser beam by an amount 
Ad, and the corresponding displacement of the focal point at the photo- 
multiplier tube is represented by At, the relationship between Ad and 


Li will depend on the focal lengths of the lenses, the diffraction 
property of the hologram, and the geometry of the system. If a linear 
relationship can be assumed and If the angle between the bisector of 
the transform lens (the broken line in Figure 3) and the front surface 
of the object is e , then the relationship between li and Ad can be 
written as 


2f c + f T 

A£ = Ad( — ^ — 4 cos(e - e 1Q ) , 


(26) 


where f c and fy are the focal length of the correlation lens and the 
transform lens respectively 





IV. Conclusions 

A theoretical analysis has been presented for the correlation 
output signal for the MSFC holographic correlation filtering system. 

It was found that under appropriate assumptions, the correlation output 
may be written as a function of the roughness of the test surface, the 
displacement of the illuminated area of the test surface, the charac- 
teristics of the optical components involved in the system, and the 
system configuration. 

Although no detailed comparison between the experimental result 
and the theoretical prediction has been made at the conclusion of this 
study, a few observations can be madei First, the theoretical result 
predicts that the correlation function has a Gaussian -form variation 
for the in-plane displacei^ent and hence no secondary peak should appear. 
The present experimental result found no secondary peak whatsoever. 
Second, the holographic filter served in a statistical averaging manner, 
and hence its exact replacement should have no significant influence 
on the output signal, this agreed with the preliminary experimental 
results as well. Third, if an aperture of a size smaller than the 
object beam crosssection is inserted in the object beam of the system, 
it can be visualized from the theoretical analysis that its effect is 
to reduce the intensity of the output signal, since the total amount of 
light diffracted by the hologram is correspondingly reduced by the 
aperture . 

Finally, if the variations of the parameters in the system are 
not sufficient to cause a variation of the output, an approximate rela- 
tionship between the displacement of the detected signal (a focused spot) 
and that of the illuminated spot on the test surface was also derived. 


Experimental verifi cation of the predicted linear relationship should 
be made in future studies on this system. 


APPENDIX A 

Deflation of Equation (14) 


If the second moment of the roughness is given by 


<ti ?2> a ° 2 C (xi» Yi» x 2 - \ (Al) 
C being the autocorrelation coefficient of the roughness, we can apply 
the lemma 


<H’(f)K'(f + g)> - -d 2 <H(f)H(f + g))/dg 2 , 


and obtain 

<^xi 5x2) = 

and 


= _ n 2 e^C 
6i x 


<cyl ty2 > = '° 2 % ' 


A x = X 2 - Xi - 

Ay = y 2 - yi - 2^ 


(A2) 

(A3) 


(A4) 


Separating the ^-dependent part of the term in angular brackets in 
Eq. (12), we then find by standard calculation the characteristic 
function 


X (v 2 i, - v z2 ) = <exp[i(v zl ? x - v z2 ? 2 )]> 

s exp(-a 2 Av z /2)e" z ^’ c \ (A5) 


where we have written Av z = v z2 - Vzi and 

g * o 2 v 21 v Z2 . (A6) 

For ordinary rough surfaces the rms roughness a is of the order 
of the wavelength x or larger, and there holds g >> 1 (in fact, g - 100). 
Substantial contributions to the integral Eq. (12) can then come only from 
those aperture elements for which C 1. We assume that C has the form 
C = exp [-( a 2 + aJ)/X 2 ] , (A7) 

where X is the correlation distance of the roughness, A x and Ay have been 
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given in Eqs. (A3) and (A4). The function 

(A8) 

with C given by Eq. (A7) can be expanded about the saddle point coor- 
dinates^ Ax * Ay * 0 or 

(A9) 

in the plane formed by 


Es » \ 


ns “-^ 


C = x 2 - Xj, n = y 2 - yi . (A10) 

With 

q 2 = g/X 2 (All) 

one obtains for the first two terms 


p ■ 1 - q 2 (a 2 + Aj) + ••• , (A12) 

which will be approximated by 1 ' 

p * exp[-q 2 (4 + A|)] . (A13) 

For very large values of g this approximation is always very good. 

For 2) » r 0 the surface elements for which C * 1, are outside 
the illuminated area, and contributions to C»i do not arise because of 
the pupil function P(x,y), Eq. (7). 

To calculate the average 
<F,F 2 >= <F(Cxi. tyi)F(cx 2 » Cy 2 ),^ 

in Eq. (12), we follow an approach already applied by Hagfors to the 
backscatter case. In essence, we assume that most of the surface slopes 
of importance with a particular scattering vector have values close to 
the mean values e x and ly defined c x = -v x /vz and ly = -Vy/vz. This 
means that Fx F 2 has a maximum for c x i = 5x2 = 5x> 5yi = 5y2 “ 5y* 
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Expanding F 1 F 2 about Cx and cy and substituting Ol in the second moments, 
Eqs. (A3) and (A4), we find that all terms but the zero-order term [F(Ex* 
Cy)] 2 are very small and can be neglected. In the calculation of the 
average we can replace e lt e 2 , 63 by the center values 6 i 0 , e 20 , 030 * 
respectively. We find 

<FiF 2 )=<F 2 ) “ 4\h 2 )(l + cos0 lo cose 2o 
- sine lo sin0 2 ocos03o) 2 / 

(cos 6 io+ cose 2 o ) 2 . (A14) 

Here R 2 is the mean square of the reflection coefficient. For the back- 
scatter case (0 2 tr-0io» 030 = 0 ) one obtains the cos" 2 0 lo dependence found 
by Hagfors? 

Substituting Eqs. (A3)-(A5), (A8), (A13), and (A14) into Eq. (12) 
and setting 

JL- k cos8 10 ((F 2 > )fy(2ir 2 r£R 0 ) , (A15) 

we can obtain the result of Equation (14). 
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